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ABSTRACT 



With IRAM-30m/HERA, we have detected CO(2-l) gas complexes within 30arcsec (~ 100 pc) from the center of M 31 amounting to 
a minimum total mass of 4.2 x 10 4 M Q (one third of the positions are detected). Averaging the whole HERA field, we have shown that 
there is no additional undetected diffuse component. We show that the above gas detection is associated with gas lying on the far side 
as no extinction is observed in the optical, but some emission is present on infra-red Spitzer maps. The kinematics is complex. (1) 
The velocity pattern is mainly redshifted: the dynamical center of the gas differs from the black hole position and the maximum of 
optical emission, and only the red-shifted side is seen in our data. (2) Several velocity components are detected in some lines of sight. 
Our interpretation is supported by the reanalysis of the effect of dust on a complete planetary nebula sample. Two dust components 
are detected with respective position angles of 37 deg and -66 deg. This is compatible with the superposition of the (PA=37 deg) disk 
dominated by the 10-kpc ring and the inner 0.7-kpc ring detected in infrared data, which position angle (-66 deg) can be measured for 
the first time. The large scale disk, which dominates the HI data, is very inclined (i=77 deg), warped and superposed on the line of 
sight on the less inclined inner ring. The detected CO emission might come from both components. 

Key words. Galaxies: individual,Galaxies: nuclei, Galaxies: kinematics and dynamics, Galaxies: bulges, Galaxies: ISM, ISM: plan- 
etary nebulae: general 



1. Introduction 



is usually described 
little star formation, 



as a quiescent galaxy 
at a level of O.4M yr _1 



M31 
with 

(e.g. Barmby et al. 2006; Tabatabaei & Berkhuiisen 2010; Azimlu 
and with ultra- weak nuclear activity ddel Burgo et alj 2000) 
The presenc e of a very massive black hole 
(Dressier 1984) and the lack o f gas with in 300 pc 
dNieten et alj|2006l: IChemin et al.ll2009t iBraun et al.l 2009b 
suggest that the main gas reservoir has been accreted and is 
exhaus ted, although some gas is detected within 1 kp c from the 
center dMelchior et al.ll2000t [M elchior & Combesl 20 1 II) From 
optical emission lines iJacobv et al.l ( 1985) estimated an ionized 
gas mass on the order of 1500M o , which can be accounted 
for by mass loss from evolving stars. iGroves et al.l (2012) rely 
on Herschel data to argue that the dust properties are well ac- 
counted for by the stellar heating. Small amounts of molecular 
gas have been detected in directions more than 300 pc from 
the c enter. These can be associated with dust featur es in this 
area dMelchior et al.ll2000l iMelchior & Com bes 201 1 ) While 



the center of M3 1 hosts a supermassive black hole with a mas s 
of 0.7 - 1.4 x 10 8 M o dBacon et al.ll2001l;lBender et all 2 0051) 
it is one of the most silent ones dGarcia e t al. 2010 ) , al - 
r^rvr "t ? mnm 8 m 2008 it started to murmur dLietal.1 20111) 

,Xa .At U. J_ i- X _ avliiV\ttc m o t-il 7 /irvliamnt f'tni/»tiirat' i f\ torn i~<^\ ta /\ 



Table 1. Log of observations 



Date 


< T w > 


^integration 


# pixels (scans) 


8 Nov. 2011 


320K 


276min 


36 (144) 


10 Nov. 2011 


411K 


276min 


36 (72) 


27 Nov. 2011 


274K 


276min 


18 (36) 


12 Feb. 2012 


25 IK 


144min 


36 (144) 


24 Feb. 2012 


318K 


300min 


36 (72) 


11 Mar. 2012 


351K 


295min 


36 (72) 



th 

e, it exhibits many coherent structures interpreted 
as tra cers of its me rging history: there is a lopsided nuclear 
disk (iLauer et a l. 1993) with two stellar components, PI and 
P2 separated by 0.45" in the center. From the kinematics, 
the black hole is located in betwe en PI and P2, but closer to 
P2. An A-star cluster (see also lKormendv & Ben der 1999) 
det ected in a third com ponent (P3) of M31's double nucleus 
bv iBender et al . (2005), can be associated with a recent star 
formation episode. This occurred 200 Myr ago, involved a total 
mass in the range 10 4 - 10 6 M o , and corresponds to an accretion 
rate of 10~ 4 - 10 _2 M Q yr _1 . Its presence so close to the black 
hole raises a number of issues: how young stars were formed 
deeply inside to the tidal field of a supermassive black hole and 
how such stars have formed while there is no cold gas detecte d 
in the surroundings (e.g. lLauer et al ] |20 1 2t iLietaD 20091) ? 
In the Galaxy, SgrA* has experienced X-ray flares, attributed 
to the infall of gas , while a cloud of gas identified by 
iGillessen et al. (2012) is expected to fall onto the black hole 
in 2013 . M31* is experiencing a similar murmur according to 
iLiet al.l (2011h suggesting some gas infall. 

Beside the presence of a young star cluster, an ionized 
gas outflow is detected in X-rays alon g the minor axis of the 
galaxy by iBogdan & Gilfanovl (2008i). perpendicular to the 
main disc. The relative intensity of the outflow on both sides is 
compatible with the intensity of the observed B extinction: the 
NW side is more extinguished than the SE side. As discussed in 
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Fig. 1. Field M31-la centered on RA: 00:42:44.1 and DEC: +41:15:42 (J2000) observed with IRAM-30m/HERA. The center is 
indicated with a cross and the circle displays the 12" beam. The Y scale of each spectrum is in main-beam temperature (T m /, 
between -0.01 and 0.02 mK, as indicated on the top right corner), while the X-axis shows velocity for the range between -600 
and kms^'and smoothed to 13 km s _1 velocity resolution. A thin line indicates the systemic velocity at -310 kms -1 . The color 
coding (displayed on the top right corner) of the velocity is used for the spectra with 3<x detections. 



iMelchior & Co mbes (2011), the velocity field of the circumnu- 
clear region (40"x40" or 150pcx 150pc) measured in optical 
ionized gas does not exhibit any clear rotation pattern: beside a 
spot at the systemic velocity in M31's center, the whole area is 
blueshifted with respect to the systemic velocity. This coherent 
flow of ionized gas is decoupled from t he stellar kinem atics 
dBender et al. 2005; del Burgo et alfeOOOt ISaglia e"t~aT]|20 101)1 
and could be connected to the recent star formation activity. 
Please note that we assume throughout the paper a distance 
to M31 of 780kpc dVilardell et all 20061) i.e. 1 arcsec=3.8pc. 



Most up-to-date r esults, based on cepheids, quote 752+27 kpc 
(Riess et al. 2012) However, for coherence with previous works 
we keep 780 kpc, which lies within l<x uncertainties of these 
new results. 

Some cold gas is expected to feed the black hole, even 
though in contrast with the Milky Way, a general lack of HI in 
the vicinity of M3 1 's nucleus has been noted for several decades 
(lEmersonlll976HBaiaia & Sha ne 1982) In this paper, we present 
the first molecular detections within 30 arcsec from the center. In 
Section|2] we present the new observations performed at IRAM- 
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Fig. 2. Field M31-la centered on R A: 00:42:44.1 and DEC : +41:15:42 (J2000) observed at IRAM- 30m HERA, superim posed on 
the A^i . observed extinction map (Melchior et al. 2000) (left panel), He and [Nil] emission map (ICiardullo et al.l 1988) (middle 
panel), and 8/imSpitzer emission map (Blo ck et al.l 2006) (right panel). The circles correspond to the > 3<x detections discussed in 
this paper. 



30m near the center of M3 1 . In Section 13.11 we analyze our 
molecular gas detections and compare them with other wave- 
lengths. In Section|U we discuss how the interpretation of these 
data.. 



2. Observations 

In the period between November 2011 - March 2012, we used 
the 1.3mm multibeam HE terodyne Receiver A rray (HERA) at 
the IRAM-30m telescope dSchusteret al.ll2004l) to do a CO(2- 
1) survey of M31's 0.7 kpc inner ring. One of the fields of 
this survey contained the center of M3 1 . We refer to this field 
as M31-la and show a 60"x60" map with 12" angular reso- 
lution in Figures Q] and [2] We thus performed a 60" by 60" 
map with 12" spatial resolution for the CO(2-l) line. Data 
were acquired in wobbler switching mode, using the Wideband 
Line Multiple Autocorrelator (WILMA) facility as backend. The 
wobbler throw (ranging between 60" to 210" in azimuth) has 
been adjusted every 30 min, to avoid extinction areas (and possi- 
ble associated molecular gas emission). This procedure has been 
followed as best as possible given the constraints of the Pool 
observing periods allocated to this project. Some signal in the 
OFF has been detected in some scans, but as the OFF positions 
are changing with time, on average it disappears in the reduc- 
tion process. Some signal can be underestimated but this should 
be included in the global 20% calibration uncertainties. A to- 
tal of 540 raw spectra were recorded in 21.5 h of telescope time 
with a spectral resolution of 2.6 kms _I . Table Q] displays the 
epoch of observations, for each day the average system temper- 
ature, the total integration time accounting for two independent 
polarization measurements and the number of pixels and scans 
performed. We provide main-beam temperatures throughout this 
paper with beam and forward efficiencies: B c yy =63 and F e //=94. 

Each spectrum was first visually inspected to remove those 
affected by very unstable baselines. They are then processed in 
order to iteratively correct for the instrumental response: (1) a 
linear baseline is subtracted for each scan, (2) the scans are aver- 
aged for each position and smoothed. (3) For the spectra where 
there is a clear detection, a higher order polynomial is subtracted 
in order to optimize the baseline subtraction. The typical rms in 
each final pixel was at lcr T mb = 3.9 mK in 13 kms'channels. 
13 of the 36 positions exhibit a detection above 3<x (as displayed 
in Figure [TJ. 



Upper limit on the continuum. We have estimated the value 
of the continuum at 1 mm in the field M31-la. 4% of the spectra 
have a mean value outside the range [-0.1,0.1] mK and have been 
removed. The remaining spectra are averaged and smoothed over 
the whole bandwidth (936 Mhz). We find a continuum level of 
0.041+0.043 mK. We thus derive a 3cr upper limit of 0. 1 3 mK on 
the continuum emission, corresponding to S cont i nuum < 0.65 mJy. 

3. Analysis 

3.1. Molecular data 

Table [2] summarizes the characteristics of the CO(2-l) 
lines detected in this field. A Gaussian function is fit- 
ted to each line to determine its area Ico, central veloc- 
ity Vb> width cr and peak temperature T pec ^. The base- 
line rms is provided for each line. We assume a standard 
Gal actic X r o =N h JI C q =2.3xl0 20 cm- 2 (K kms" 1 )- 1 follow- 
ing IStrong et al.l ( 1988) However, no te that diffe rent values 
have bee n adopted in the literature. iLero v e t aD ( 20 1 1 1) re- 
lying on iNieten et al.l (20061) CO data with a strong signal 
(Jco >lKkms _I ), estimate a lower value Xqo = 9.66 + 1.33 x 
10~ 19 cirr 2 (,K km s~ l )~ l for the inner part of M3 1 . 

Relying on the CO(2-1)/CO(1-0) line ratios measured by 
iMelchior & Co mbes (2011) in this area, we assume a line ratio 
of 1 and thus adopt the previous Xco ratio for the CO(2-l) line. 
We then convert the N#, column density to an H2 mass surface 
density and derive a molecular mass M/ 3f , flm assuming the gas 
fills the main-beam. Lastly, when all the positions are averaged 
a noise level of 0.8 mK is achieved, but no signal appears. This 
stacking demonstrates that there is no extended emission larger 
than 2.4 mK. 

One third of the observed positions exhibit a CO(2-l) detec- 
tion, as displayed in FigureQ] Beside the spectra 18 and 23 (and 
the second component of spectra 36), all the detected lines are 
redshifted with respect to the systemic velocity. 

3.2. Dust extinction and dust emission 

In this area devoid of large amounts of gas 
(e.g. Loinard et al. 1996), we have shown in 
Melchior & Combes (2011) that in the North- Western part 
of the bulge of M31, CO is detected where extinction is 
observed, while it is not detected in areas where no extinction 
is measured. This supports the dust-gas correlation observed 
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Table 2. Characteristics of the CO(2-l) lines. Spectra are smoothed to a 13 kms 'resolution, but the fits in italic have been per- 
formed on 2.6 km s spectra. The offsets refer to the center of the M3 1-1 a field. 



# 


Offsets 


R(") 


Ico (Klaus -1 ) 
= \T mb dV 


V (km s -1 ) 


cr (km s ') 


rppeak 

mb 

(mK) 


baseline 

rms 

(mK) 


N H2 (cnT 2 ) 


(M G pc~ : 


M beam (M Q ) 

) 


5 


-33.0,-8.33 


34.0 


0.51+0.12 


-301.7+4.0 


30.1+8.6 


16.0 


4.2 


1.17 x 10 20 


1.99 


3.25 x 10 3 


15 


-9.0,-32.3 


33.5 


0.60±0.08 


-248.3+1.7 


24.4+3.9 


23.0 


3.2 


1.38 x 10 20 


2.35 


3.84 x 10 3 


17 


-9.0,-8.3 


12.2 


0.76±0.13 


-137.4+3.9 


41.9+7.3 


17.1 


4.0 


1,75 x 10 20 


2.97 


4.85 x 10 3 


17 


-9.0,-8.3 


12.2 


0.26±0.08 


-32.7+3.9 


13.3+29.1 


18.4 


4.0 


0.60 x 10 20 


1.02 


1.67 x 10 3 


18 


-9.0,-3.7 


9.7 


0.35±0.09 


-355.0+4.1 


27.8+7.0 


11.7 


3.5 


0.81 x 10 20 


1.37 


2.24 x 10 3 


19 


-3.0,-56.3 


56.4 


0.13+0.05 


-341.4+4.0 


13.1+3.9 


9.2 


3.8 


0.30 x 10 20 


0.51 


0.83 x 10 3 


19 


-3.0,-56.3 


56.4 


0.10±0.03 


-339.9+4.0 


7.5+2.2 


12.9 


5.1 


0.23 x 10 20 


0.39 


0.64 x 10 3 


21 


3.0,-32.3 


32.4 


0.50+0.13 


-224.0+6.9 


48.1 + 11.1 


9.8 


3.9 


1.15 x 10 20 


1.96 


3.20 x 10 3 


22 


3.0,-20.3 


20.5 


0.20+0.07 


-67.1+4.0 


17.2+6.2 


11.0 


3.0 


0.46 x 10 20 


0.78 


1.27 x 10 3 


23 


3.0,-8.3 


8.8 


0.34+0.09 


-390.4+4.0 


29.5+8.1 


10.9 


3.2 


0.78 x 10 20 


1.33 


2.17 x 10 3 


24 


3.0,3.7 


4.8 


33+0 09 


-256 6+5 


28.6+8.6 


11.0 


3.1 


0.76 x 10 20 


1.29 


2.11 x 10 3 


26 


15.0,-44.3 


46.8 


0.48+0.14 


-207.8+5.7 


39.8+14.5 


11.2 


4.1 


1.10 x 10 20 


1.88 


3.07 x 10 3 


26 


15.0,-44.3 


46.8 


0.45+0.16 


-121.1+7.7 


47.0+21.9 


9.1 


4.1 


1.03 x 10 20 


1.76 


2.87 x 10 3 


28 


15.0,-20.3 


25.2 


0.36+0.08 


-190.9+1.0 


13.0+33.1 


26.1 


3.8 


0.83 x 10 20 


1.41 


2.30 x 10 3 


28 


15.0,-20.3 


25.2 


0.22±0.05 


-189.7±1.9 


14.1 ±3.1 


14.4 


5.6 


0.51 x 10 20 


0.86 


1.40 xlO 3 


35 


27.0,-8.3 


28.2 


0.15+0.05 


-160.9+3.3 


13.0+4.79 


10.9 


3.6 


0.35 x 10 20 


0.59 


0.96 x 10 3 


35 


27.0,-8.3 


28.2 


0.064+0.026 


-155.0±0.8 


3.4±1.8 


17.8 


5.4 


0.15 x id 20 


0.25 


0.41 x 10 s 


36 


27.0,3.7 


27.3 


0.94+0.16 


-152.8+4.6 


52.6+10.6 


16.8 


4.1 


2.16 x 10 20 


3.68 


6.01 x 10 3 


36 


27.0,3.7 


27.3 


0.43+0.10 


-315.6+2.8 


21.9+5.3 


18.3 


4.1 


0.99 x 10 20 


1.68 


2.74 x 10 3 



in t he Milky Way dBohlin et al . 1978) and other galaxies 
(e.g. lFovle etall 20121) ISmith et al] (20121) is also claiming 
that the gas in M31 is well traced by dust at a constant 
metallicity. The optical and near-infrared data displayed in 
Figure [2] provide complementary information to our CO 
detection. T he left panel displays the observed extinction as 
computed in Melc hior et al.l (2000) The bulge light is mostly 
dominant within R < 1.2kpc (~ 300arcsec) from the center 
dCourteau et al J 20 1 II) j where R is the projected distance to the 
center on the sky plane. We model its photometry with elliptical 
annuli using the standard surface photometry algorithm devel- 
oped for IRAF dJedrzei ewski 1987) This model is intended 
to reproduce the light profile along the bulge of M31 without 
extinction. The median intensity over the elliptical annulus 
sectors is used, in order to avoid areas suffering extinction. 
Large scale extinction following the elliptical profile could affect 
this model but it is obviously weak. The observed extinction 
is defined as A A = -2.5log l0 (<f> obs /(p mode i), where <p obs is the 
observed brightness and <p mo dei the modeled brightness. This 
provides a lower bound of the true internal extinction, with the 
following assumptions: (1) the gas lies in front of the bulge 
and (2) there is only one clump of gas per line of sight. It thus 
eliminates automatically large scale homogeneous extinction 
like the one due to the Milky Way. On the one hand, as the 
resolution of the optical data is close to 1 arcsec corresponding 
to 3.8 pc at the distance of M3 1 (the optical resolution is a factor 
of 10 better than the one achieved with CO observations), the 
second assumption is reasonable for most lines of sight given 
the small masses detected. According to the mass-size r elatio n 
for molecular clouds first observed by Solomon et al. (1987), 
we expect that the more massive molecular clouds detec ted here 
have a size smaller than 3pc (e.g. lLombardi et aUl2010l) . On the 
other hand, it is most probable that all the gas does not lie in 
front of the bulge in the Southern part. The observed extinction 
corresponds to A observed = -2.51og 10 [/ + (1 - /) x exp(-T)], 
where / is the fraction of light in front of the dust and t 
the real optical depth at a given wavelength. In principle, if 



several clumps were present along the line of sight, we could 
decompose the extinction A observa i into several components. In 
practice, more observational constraints would be necessary to 
perform such a detailed analysis. It is important to note that 
the majority of the dust and associated gas is not diluted in 
the whole bulge and not spread all along the line of sight. It 
is expected to be relaxed and to lie in a disk-like structure (or 
ring). 

The top panel of Figure[3]displays how the observed extinc- 
tion relates to the fraction / of light in front of the dust for 
different values of real internal extinction. For extinction val- 
ues above 5, it is very difficult to disentangle the real value of 
the extinction from optical observations. For a given value of 
observed extinction, it is, however, possible to put some con- 
st raint on the po sitions of dust clumps. Relying on the modeling 
of iTempel et al.l (201 II) we estimate how the fraction / of light 
in front of the dust relates to the position z of the dust clump 
along the line of sight. The bottom panels of Figure [3] display 
this relation for different lines of sight along the minor axis Y 
(projected on the sky plane). If the line of sight crosses the cen- 
ter (where there is a cusp due to the stellar cluster), the fraction 
/ is varying strongly as most of the light comes from an area 
smaller than 1 arcsec (~ 3.8 pc). For lines of sight varying from 
10 pc to lkpc, the slope of f(z) is decreasing. There are two 
main effects to be stressed. (1) The observed extinction is very 
sensitive to the location of the dust clump along the line of sight. 
In the central region, dust just behind the mid-plane could easily 
escape optical detection, as the fraction of light in front of the 
clump can be very important. (2) This effect is further strength- 
ened by the asymmetry due to the inclination: the far side will 
have a larger fraction of light than the near side. 

The real extinction is difficult to measure from optical data, 
but the observed extinction provides some constraints on the lo- 
cation of the dust along the line of sight as displayed in Figure|2] 
For instance, if Absented = 0.2, the real extinction is larger than 
0.25 and the fraction of light / in front of the dust between 0.2 
and 0.85. Then as a function of the chosen line of sight, it is pos- 
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Fig. 3. Relationship between extinction and the position of the 
dust clumps along the line of sight. The upper panel displays 
how the observed extinction A t ser ved relates to the fraction / of 
light in front of the dust for different values of real extinction 
A rea i. The lower panels display how the fraction / relates to the 
line of sight distance z with respect to the center for different 
projected distances Y to the center along the minor axis. We 
detect a near/far side asymmetry, which strengthens the effect of 
absence of observed extin ction on the far side. We rely on the 
modeling of Tempe l et al.l (2QJM 



sible to constrain the position of the dust clump. For a position at 
lOOpc (resp. lOpc) from the center, the dust clump is expected 
to lie within 400 pc (resp. 150pc) from the plane perpendicular 
to the line of sight passing through the center. Note that these 
observations have a resolution of 12 arcsec (~ 45 pc), so the gas 
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Fig. 4. Presen ce of two dust components affecting the 

ICiardullo et al. ( 1989) sample of planetary nebulae. We consider 
a complete sample of 99 planetary nebulae distributed into a 15- 
200 arcsec annulus. This annulus has been split into five equal 
parts as indicated by the colors. The left panel displays the re- 
lation between the number N of planetary nebulae detected in 
each position angle. The color star points (on the left panel) cor- 
respond to the color hatched areas displayed on the right panel. 
The black stars (on the left panel) have been obtained with a run- 
ning sum for different position angles (PA). The error bars indi- 
cated only on the color points correspond to Poisson statistical 
noise. The line corresponds to the multiplication of two sinu- 
soids of phase PA=37 deg and -66 deg, as indicated by the ticks 
at the bottom. The right panel displays the spatial distribution of 
the planetary nebula sample in the annuli corresponding to the 
completeness limits (15, 90 and 200 arcsec). 



present very close (in projected distance) to the black hole could 
lie anywhere between and 150pc on the far side. 

The right panel of Figure [2] displays the dust emission at 
8 fim. It is not affected by extinction, but depends on dust grains 
and their heating. The Southern part exhibits a dust-emission in- 
tensity much stronger than the observed extinction, suggesting 
that the dust clumps lie on the far side. 



3.3. Characteristics of the dust components 

Following the previous discussion on the expected near/far 
asymmetr y expected due to the inclination of M3 1 , one can 
note that ICiardullo et al.l ( 1989) first detected an asymmetry 
between the near and far sides among the bulge's planetary 
nebulae distribution. A similar geometrical effect is also ex- 
pecte d in the distributi on of microlensing events in M31's bulge 
(e.g Kerins et al. 2001, 2006) An asymmetry between the near 
side and the far side, due to extinction in the main plane, is ob- 
served in the distrib ution of planetary nebulae in the lOkpc ring 
dMerrett et alJ 20061) 

We have re i nvest igated the catalog of planetary nebulae of 
ICiardullo et al. ( 1989), which samples quite well the bulge area. 
There is a com plete sample of 99 planetary nebulae: it is spatially 
complete and ICiardullo et all ( 1989) computed the detection ef- 
ficiencies with respect to the surface magnitude. As displayed 
in the right panel of Figure |U 29 are present in a 15-90 arcsec 
annulus within a m5007 magnitude smaller than 22.1, and 70 in 
a 90-200 arcsec annulus with a m5007 magnitude smaller than 
22.7. This region within 200 arcsec from the center is dominated 
by the bulge, and the planetary nebulae follow the light distribu- 
tion. In the left panel of Figure [4] we have counted the number 
of planetary nebulae in five parts of the 15-200 arcsec annulus 
and displayed them as a function of the position angle. The five 
points follow a sinusoid. The star points (not independent) have 
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been obtained similarly with a running sum for intermediate po- 
sition angles. 

According to the near/far side asymmetry, one would expect 
a sinusoidal variation of the number of planetary nebulae with 
a 2n period with respect to the position angle of the main disk 
37 deg. Surprisingly, the observed period is n. The overplotted 
sinusoid is varying as sin(PA - 37) x sin(PA + 66), which is com- 
patible with the superposition of two dust components whose 
main axis have respective position angles of 37 degrees and - 
66 degrees. This is a new confirmation of the presence of two 
gas/dust components in this region: the main disc (and mainly 
the lOkpc ring) seen in projection with a PA of 37 degree, and 
the inner ring seen in the infra-red (e.g. with Spitzer data) with a 
position angle of -66 degree. The amplitude of the effect is sim- 
ilar for both components. There is probably an additional per- 
turbation (at the limit of detection) close to PA=- 180/1 80 deg, 
which might correspond to a non-circular structure. More statis- 
tics are required to be more conclusive. 

Those two dust components expected to be associated with 
CO detections have an orientation compatible with our previous 
discussion: the clumps detected in CO are most probably on the 
far-side. In addition, one could argue that the few points detected 
close to the systemic velocity (and blue shifted) could be associ- 
ated to the main disc (possibly at large scale).. 

3.4. Ionized gas 

The middle panel of Figure [2] show s the Ho- and [Nil] 
emission map (ICiardullo et al.l 19881) As discussed by 
Rubi n & Fordi d 97 11) it is dominated by [Nil] excited by 
shocks. The overall pattern corresponds to the 8 /urn map, 
but there is not an exact correspondence. (1) These wave- 
lengths are affected by extinction. For ins tance, the position 
observed in CO by iMelchioret al. (2000) (double circle in 
Figure [2]i seems affected by extinct ion. (2) The kinemat- 
ics of t he [Nil] line measured by iBou lesteix et ail ( 19871) 
(see also Melc hior & Combesl 201 1) exhibit a disc in rotation 
and the (40"x40") circumnuclear region is blueshifted with 
respect to the systemic velocity. In parallel, the velocity field 
measured in CO is redshifted. It is thus probable that both 
components are decoupled. 

This can be related to the co mpilation of all gas velo cities 
in the inner lOarcmin of M31 by IStark & Binnevl d!994t) . The 
isovelocity curves are very irregular and chaotic, even involving 
much larger scales than here. 



3.5. Bulge light emission 

The bulge model computed on photometric images as de- 
scribed in Sec tion 13. 21 provides the ellipse geometry parameters 
(see e.g. Kent 1989) We performed this modeling on 2MASS J 
datcQ (Skrutskie et al. 2006) and Figure displays the position 
of the centers of each annulus computed. In contrast with B and 
Ha/[NII] data, we do not expect any bias due to dust obscuration. 
Interestingly, the center of the annuli is systematically shifting 
towards the South by about 7pc within 0.78 kpc. It is tempting 
to compare this off-centering of the bulge w ith the off-centerin g 
(~ 350 pc) of the inner dust ring detected bv lBlocketall (20061) 
Under the hypothesis of a coupled m - 1 motion betwee n the 
inner bulge and the disk, and given the ITamm et al. (2012) esti- 
mation the mass of M31's bulge of the order of 4 x 10 I0 M o , we 
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Fig. 5. Variations of the position of the centers of the elliptical 
annuli computed on the 2MASS J image as a function of the 
semi-major axis. 



can expect a mass of gas in the inner ring of 9 x 10 8 M o , if the 
maximum of light corresponds to the barycenter. 

As discussed in the Appendix[A]and shown in Figure lATTl the 
bulge is triaxial but we do not expect our results to be affected as 
the amplitude of the twist is about 10 deg. 



4. Interpretation 

As displayed in Figure [2j some detections (7/13) correlate well 
with the Spitzer dust emission and several position^ (5, 17, 18, 
19, 23 and 26) do not. None of the positions detected in CO cor- 
respond to any observed A B , , extinction. According to the 

* J observed ° 

available kinematics, it is most probable that the ionized com- 
ponent is decoupled from the molecular gas. The observed ex- 
tinction is smaller than 0.025 in B (see left panel of Figure 0. 
Relying on Figure[3] we can derive that the typical fraction / of 
light in front of t he average clump is larger than 90%. Following 
the modeling o f lTempel et ail (20111)] the average clump lies at 
a depth between 20 and 200 pc from the center on the far side, 
depending on its projected distance. It could be further if the 
real extinction is significantly larger than 0.25. Accordingly, the 
positions with no Spitzer infrared emission could be much fur- 
ther where the radiation field is too weak to heat the dust, except 
within 4 pc of the black hole, where the light of the nuclear star 
cluster prevents the detection of any extinction. 

The kinematics is complex and do not exhibit a clear pattern. 
Most of the lines are redshifted with respect to the systemic ve- 
locity while 4 are blueshifted. This is surprising as the opposite 
trend was observed in this field bv lMelchior & Combesl (201 If) 
with the optical ionized ga s. The gas detected in the North-East 
side dMelchior et all 2000) was also redshifted, so it is not sim- 



1 Note that the J image has been chosen as it offers a better signal-to- 
noise ratio than the K image. 



2 The Spitzer map seems to exhibit a defect close to the center at the 
position of spectrum 17. 
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ply a counter-rotation. In addition, the velocity range is spread 
between -33 and -390 kms -1 . 

These results a re compatible with the 0.7 k pc inner ring sce- 
nario discussed in Melchior & Combes (201 1) and initially pro- 
posed by Bloc k et alj (20061) : the ring is tilted, which could ex- 
plain why the gas is here on the far side and the velocities red- 
shifted. It lies off-center, which could explain why we do not 
see a regular rotation pattern, since we are far from the kine- 
matical center. As supported by the dust components detected 
with the planetary nebulae in Sect. 13.31 the 0.7-kpc inner ring is 
most probably superimposed on the 10-kpc ring in the main disc: 
this could account for the multiple velocity components together 
with the dumpiness. 

In summary, we have shown the presence of molecular gas 
close to the black hole. There is no extended diffuse molecu- 
lar emission but we have detected small dense clumps located 
on the far side of the bulge. The detected clumps are located 
between 20 pc and 215 pc in projected distance from the cen- 
ter (and observed with a resolution of 45 pc). Assuming a single 
dust/gas clump per line of sight and some modeling assump- 
tions, we show that clumps corresponding to these lines of sight 
lie on the far side, at least at 150pc in depth from the center and 
most probably closer than 600 pc. If some gas is present next 
to the center, it is also most probably on the far side but could 
be very close to the black hole. The kinematics exhibit an un- 
expected trend: most of the lines are redshifted with respect to 
the systemic velocity, which could be due to the off- centering. 
Along several lines of sight, there ar e m ultiple CO components. 
Our reanalysis of the Ciardullo et al. ( 1989) catalog of planetary 
nebulae reveals the presence of two components of position an- 
gles 37 deg and -66 deg. In the light of our molecular detections 
discussed here, we claim that dust and gas from the outer 10-kpc 
ring superimpose on the inner 0.7-kpc ring in this very central 
region. The inclination of both components is such that the gas 
and dust in the Southern area lie mainly on the far side. 
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Appendix A: Triaxality of the bulge 



Figure lATTl displavs the variation of the position angles and ellip- 
ticities computed in the modeling described in Sect. 13.21 There 
is a clear isophot twist which is not due to extinction. The most 
pl ausible expla nation is the triaxiality of the bulge as discussed 
bv lKormendv (e.g. 1982 5, Ou r values are similar to those pre- 
sented bv lBeaton et alj (20071) , but they are presented here in lin- 
ear scale in accordance with Figure [5] 

The isophot twist is significant but its amplitude does not ex- 
ceed lOdeg and should not affect significantly our results. (1) 
The bins used to compute the position angle of the dust compo- 
nent from the planetary nebulae distribution are 72 deg. (2) This 
triaxiality could affect the near/far side effect, but it should be a 
second order effect. 
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Fig.A.l. Position angles and ellipticities of the centers of the 
elliptical annuli computed on the 2MASS J image as a function 
of the semi-major axis. These points correspond to the centers of 
the annuli presented in Figure [5] 
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